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摘要 :【 目的] 驱 避 剂 是 防止 蚊虫 叮咬 的 良好 个 人 防护 用 品 ,但 其 作用 机 理 一 直 未 被 解释 清楚 。 本 研究 利用 化 学 计 
算 手 段 对 4 种 重要 的 人 体 引诱 物 ( 乙酸 丙 酸 、 丁 酸 和 戊 酸 ) 与 22 个 情 类 驱 避 化 合 物 的 缔 合 作用 以 及 该 作用 对 驱 避 
活性 的 影响 进行 了 研究 ,以 期 为 蚊虫 驱 避 机 理 研究 提供 新 思路 。【 方 法 】 利 用 GAUSS View 软件 构建 引诱 化 合 物 分 
子 、 驱 避 化 合 物 分 子 以 及 二 者 可 能 形成 的 缔 合 体 结构 ,再 利用 GAUSSIAN 软件 优化 所 有 结构 ,并 计算 它们 的 能 量 和 
各 类 参数 ;借助 AMPAC 和 CODESSA 软件 计算 出 所 有 结构 的 描述 符 ; 以 莫 类 驱 避 化 合 物 对 白 纹 伊 蚊 4edes albopictus 
校正 驱 避 率 的 对 数值 (logCRR ) 为 活性 数据 ,利用 CODESSA 软件 建立 描述 符 与 驱 避 活性 的 定量 构 效 关系 模型 。 
【结果 】4 种 羚 酸 与 音 类 驱 避 化 合 物 相互 缔 合 的 缔 合 能 量 大 多 在 18 ~20 kJ/mol 之 间 , 缔 合 距离 和 角度 分 别 在 2.0 ~ 
2.5 A 和 110° ~ 180* 内 波动 ,符合 氢 键 相互 作用 的 特征 。4 个 最 优 定量 构 效 关系 模型 的 尼 本 大 0.9 ,模型 中 的 
参数 都 包含 缔 合 体 的 描述 符 ,表明 引诱 物 - 驱 避 物 缔 合 作用 【结论 ]4 种 羚 酸 与 音 
类 驱 避 化 合 物 之 间 存 在 中 等 强度 的 氨 键 缔 合 作用 ,该 缔 合 作用 对 驱 避 活性 具有 重要 影响 。 0 二 未 下 
物 - 驱 避 物 缔 合 作用 的 存在 ,也 为 从 新 的 角度 研究 蚊虫 驱 避 机 理 提 供 了 前 | 
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Abstract:【 Aim】Mosquito repellents are a group of compounds which can protect human from mosquito 
biting. However, the repelling mechanisms are still unclear. In this study ，theoretical calculations were 
applied to investigate the molecular association between four important human-derived attractants (acetic 
acid, propionic acid, butyric acid and pentanoic acid) and twenty-two terpenoid repellents, to find out 
the relationship between this kind of interaction and repellency, and so to clarify the repelling 
mechanisms from a new perspective. [ Methods 】 Structures of attractants, repellents and attractant- 
repellent complexes were built and optimized by GAUSS View and GAUSSIAN, and association energy 
and geometrical parameters were obtained. Descriptors for model building were calculated by AMPAC and 
CODESSA. The quantitative structure-activity relationship ( QSAR) models were built to analyze the 
relationship between the structural descriptors and the logarithm value of the corrected repelling time 
(logCRR) against Aedes albopictus. [Results) The association energy between four carboxylic acids and 


twenty-two terpenoid repellents is between 18 and 20 kJ/mol, and the association distance and angle 
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are fluctuating from 2. 0 to 2.5 A and from 110° to 180°, respectively, which exactly fit the 
characteristics of hydrogen bonding interaction. Four obtained QSAR models were statistically satisfying ， 
and the R’ values of any one of these four parametric models were above 0.9. Every model contained the 
descriptors generated from attractant-repellent complex, showing a significant correlation between 
attractant-repellent association and the logCRR. 【 Conclusion 】 Hydrogen bonding associations with 
moderate strength exist between four carboxylic acids and terpenoid repellents, and such kinds of 
associations deeply affect the repellency. This study provides a preliminary support to the existence of 
attractant-repellent association and sets up a foundation for the research of repelling mechanism from a 
new perspective. 


Key words: Terpenoid; mosquito repellent; carboxylic acid; molecular association; repellency; 


quantitative structure-activity relationship (QSAR ) 


驱 避 剂 所 含 的 驱 避 化 合 物 是 一 类 由 植物 产生 或 
人 工 合 成 的 具有 驱 避 昆虫 作用 的 活性 化 学 物质 。 蚊 
虫 的 叮咬 不 仅 使 人 厌烦 ,同时 还 会 传播 多 种 严重 的 
疾病 。 蚊 虫 驱 避 剂 的 使 用 由 于 可 以 有 效 减少 蚊虫 叮 
咬 , 从 而 在 减少 蚊 媒 疾病 方面 具有 重要 作用 ,使 得 驱 
避 剂 的 开发 凸显 重要 。 随 着 驱 避 剂 开 发 工作 的 深 
入 ,更 为 全 面 地 了 解 驱 避 作 用 机 理 对 于 有 目的 性 地 
开发 高 效 驱 避 剂 显得 尤为 关键 。 

经 过 数 十 年 的 研究 , 驱 避 作用 机 理 研究 已 有 了 
较 大 进展 。Davis 和 Rebert(1972 ) 发 现 蚊虫 触角 上 









































乎 表明 还 有 其 他 补偿 机 制 存在 ;其 三 ,目前 的 研究 主 
要 集中 于 蚊虫 的 触角 部 位 ,对 驱 避 剂 使 用 部 位 的 研 
究 不 足 。 男 外 ,Dogan 等 (1999 ) 还 通过 行为 实验 得 
出 了 DEET 对 埃及 伊 蚊 表现 出 引诱 活性 ,而 上 只 有 在 
乳酸 同时 存在 时 , 才 表 现 出 对 乳酸 引诱 活性 的 削弱 
作用 。 这 一 结果 显示 ,在 驱 避 过 程 中 , 驱 避 物 与 引诱 
物 相互 作用 是 一 个 非常 重要 的 因素 。 

关于 引诱 物 与 驱 避 物 相互 作用 (也 称 作 缔 合 作 
用 ) 的 研究 , 目前 已 取得 一 些 进展 ( 廖 圣 良 等 , 2011， 
2012; Song ef al., 2013)。 涉 及 的 人 体 引 诱 物 包括 





























存在 对 避 蚊 腕 (DEET) 敏感 的 嗅觉 受 体 神经 元 。 之 
后 ,Davis 和 Sokolove(1976 ) 提出 驱 避 剂 DEET 可 以 
抑制 位 于 埃及 伊 蚊 4edes egypii 触角 上 对 乳酸 敏感 
的 神经 元 。Ditzen 等 (2008 ) 通过 行为 和 电 生 理 实 验 
表明 , 驱 避 剂 DEET 将 冈比亚 按 蚊 Anopheles gambiae 
的 嗅觉 受 体 作 为 分 子 裔 标 ,并 且 抑 制 了 由 气味 引起 
的 嗅觉 受 体 OR-OR83b 复合 体 的 激活 反应 。 与 此 同 
时 ,Syed 和 Leal (2008 ) 发 现在 致 倦 库 蚊 Culex 
guinguefasciatus 的 触角 上 存在 对 DEET 敏感 的 嗅觉 
神经 元 ,蚊虫 在 感知 DEET 后 ,可 以 做 出 主动 躲避 行 
为 。 以 上 驱 避 作用 机 理 研究 大 致 可 以 总 结 为 两 类 模 
式 ,一 类 为 蚊虫 感知 驱 避 剂 后 主动 和 虞 避 , 另 一 类 是 驱 
避 剂 掩盖 蚊虫 对 引诱 物 的 感知 。 不 过 在 驱 避 剂 的 使 
用 过 程 中 ,还 涉及 人 体 引 诱 物 的 存在 、 驱 避 剂 分 散 、 
引诱 物 与 驱 避 剂 共存 等 多 个 环节 ,因此 驱 避 作用 机 
理 异 常 复杂 , 目前 的 研究 结果 尚 不 能 全 面 阐释 驱 避 
剂 的 作用 机 理 。 

在 以 上 机 理 获得 广泛 关注 的 同时 ,还 存在 着 其 
他 有 待 解决 的 疑问 。 首 先 ,上 述 理 论 尚 不 能 完全 人 解 
释 在 使 用 驱 避 剂 之 后 ,其 他 人 体 主 要 引诱 物 (如 :小 
分 子 羧 酸 ) 无 法 被 蚊虫 感知 的 原因 ; 其次, Stanczyk 


















































氨 、 二 氧化 碳 、 乳 酸 等 。 不 过 ,在 人 体 分 泌 的 众多 引 
诱 物 中 ,还 存在 一 系列 小 分 子 羧 酸 , 它 们 对 于 引诱 蚊 
虫 也 起 到 了 重要 作用 。Carlson 等 (1973 ) 报道 过 C,- 
C5; 的 羧 酸 与 CO 共存 时 对 埃及 伊 蚊 具有 协同 引诱 作 
用 。Bosch 等 (2000 ) 研究 表明 一 定 浓度 的 C1-C;， 
Cs-Cs 和 Cis-Cis 的 饱和 脂肪 酸 可 以 显著 增强 乳酸 对 
埃及 伊 蚊 的 引诱 力 ,并 且 低 浓度 下 的 C1-C 饱和 脂 
肪 酸 表现 出 更 为 出 色 的 增强 引诱 效果 。 因 此 ,深入 
研究 这 些小 分 子 羧 酸 与 驱 避 化 合 物 的 相互 作用 ( 缔 
合作 用 ) 对 于 驱 避 作用 机 理 的 揭示 有 着 积极 的 意 
义 。 本 人 研究 选用 了 4 个 小 分 子 送 酸 (乙酸 、 丙 酸 、 丁 
酸 和 戊 酸 ) 入 类 驱 避 化 合 物 为 研究 对 象 ,参照 已 
有 的 研究 方法 ( 廖 圣 良 等 ，2012; Song et al.， 
2013 ) , 即 利用 理论 计算 对 4 种 羧 酸 与 昔 类 驱 避 化 合 
物 的 缔 合 作用 进行 计算 ,并 利用 定量 构 效 关系 研究 
来 说 明 这 种 相互 作用 对 驱 避 活性 的 影响 。 

















1 材料 和 方法 


1.1 化 合 物 及 其 活性 数据 
本 次 研究 的 对 象 包括 4 个 小 分 子 羧 酸 ( 乙 酸 、 











等 (2010) 曾 通过 实验 表明 即使 在 DEET 存在 的 情 
况 下 ,仍然 有 一 小 部 分 蚊虫 会 被 人 体 气味 所 引诱 , 似 











丙 酸 、 丁 酸 和 戊 酸 ) 及 22 个 昔 类 驱 避 化 合 物 。22 个 
茧 类 驱 避 化 合 物 均 由 天 然 精 油 松节油 经 化 学 衍生 化 
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制备 而 得 ,这 些 化 合 物 的 合成 和 驱 蚊 活 性 测试 参见 
王 宗 德 等 (2007) 和 Wang 等 ( 2008)。 表 1 列 出 了 
这 22 个 莫 类 化 合 物 的 结构 和 了 驱 蚊 活性 (校正 驱 避 
率 的 对 数值 logCRR ) 信息 。 
1.2 法 酸 与 若 类 驱 避 化 合 物 缔 合 的 能 量 、 距 离 和 角 
度 的 计算 

在 GAUSS View 4 软件 中 ,分 别 构 建 引诱 化 合 物 
分 子 、 草 类 驱 避 化 合 物 分 子 以 及 引诱 化 合 物 与 驱 避 
化 合 物 缔 合 作用 时 (引诱 物 - 驱 避 物 缔 合 体 ) 的 三 维 
结构 式 。 在 构建 缮 合体 的 三 维 结构 式 时 ,所 有 可 能 
参与 相互 作用 的 基 团 或 片段 都 将 在 构建 过 程 中 加 以 
考虑 ,包括 :引诱 化 合 物 中 的 羟基 . 痰 基 , 驱 避 化 合 物 
中 的 普 基 羟基 和 醚 键 。 缔 合体 的 构建 过 程 通常 是 
通过 将 引诱 化 合 物 围绕 驱 避 化 合 物 立 体 旋 转 , 以 寻 
找 所 有 可 能 的 构象 。 随 后 ,在 GAUSSIAN 03W 软件 
中 完成 对 这 些 结构 的 优化 计算 ( Frisch et al.， 
2004) ,计算 水 平 为 HFA6-31G(d) ,优化 最 终 得 到 的 
是 能 量 最 低 的 优势 构象 ,同时 还 将 获得 缔 合 体 缔 合 
部 位 的 各 项 几何 参数 。 从 GAUSSIAN 03W 的 输出 
文件 中 可 以 读 取 各 个 结构 的 能 量 ,以 及 缔 合 体 中 包 
括 缔 合 距离 . 缔 合 角度 在 内 的 多 种 几何 参数 。 

引诱 化 合 物 与 驱 避 化 合 物 缔 合 作用 的 能 量 可 以 
用 以 下 公式 计算 获得 : 

E(association) = E(A) + E(R) - E(AR) 

其 中 ,E(association) 表示 缔 合 作用 能 量 ,5(A) 
是 引诱 化 合 物 的 能 量 ,A(R ) 为 驱 避 化 合 物 的 能 量 ， 
E( AR) 代 表 缔 合体 的 能 量 。 
1.3 描述 符 的 计算 

以 图 形 的 形式 来 表示 化 合 物 的 结构 是 研究 人 员 
常用 的 方式 之 一 ,然而 在 研究 物质 的 化 学 结构 与 其 
活性 的 相互 关系 时 ,人 们 通常 需要 对 化 学 结构 进行 
数字 化 处 理 , 即 利 用 结构 描述 符 来 描述 整个 结构 。 
描述 符 通常 包括 两 类 :分 子 描述 符 和 局 部 描述 符 。 
分 子 描述 符 是 利用 软件 对 整个 分 子 进行 计算 获得 的 
描述 符 , 具 体 的 描述 符 类 型 包括 6 类 :结构 组 成 描述 
符 ,拓扑 描述 符 、 几 何 描 述 符 、 静 电 描 述 符 、 量 子 化 学 
描述 符 和 热力 学 描述 符 。 局 部 描述 符 是 就 特定 分 子 
片段 或 结构 区 域 进行 计算 获得 的 描述 符 , 图 1 所 示 
的 是 4 种 羧 酸 与 1 号 化 合 物 缔 合 作用 时 缔 合 区 域 的 
选择 (红色 虚线 框 内 ) , 缔 合 区 域内 基 团 的 选取 由 缔 
合 距 离 和 角度 决定 。 局 部 描述 符 类 型 包括 :结构 组 
成 描述 符 ,拓扑 描述 符 、 几 何 描 述 符 和 量子 化 学 描述 
符 。 所 有 描述 符 的 计算 是 在 AMPAC 和 CODESSA 
中 完成 。 
























































图 1 引诱 化 合 物 - 驱 避 化 合 物 缔 合 体 及 缔 合 区 域 
Fig.1 Attractant-repellent complexes and association areas 
A: 乙酸 - 驱 避 化 合 物 缔 合 体 Acetic acid-repellent complex B: 丙 酸 - 驱 
避 化 合 物 缔 合 体 Propionic acid-repellent complex; C: 丁 酸 - 驱 避 化 合 
物 缔 合 体 Butyric acid-repellent complex; D: 成 酸 - 驱 避 化 合 物 缔 合 体 
Pentanoic acid-repellent complex. 红色 虚线 框 表 示 4 种 凑 酸 与 1 号 化 
合 物 缔 合 作用 时 的 缔 合 区 域 。The fractions in the red boxes are the 


areas of molecular associations between four carboxylic acids and number 





one compound. 


1.4 模型 的 构建 和 检验 

利用 CODESSA 软件 中 的 启发 式 回归 方法 ,在 
已 获得 的 描述 符 中 筛选 出 与 活性 密切 相关 的 描述 符 
用 于 构建 多 元 线性 QSAR 模型 。 模 型 构建 和 检验 的 
步 又 如 下 : 

1) 将 22 个 化 合 物 的 活性 值 (logCRR ) 按 降 序 
排列 ; 

2) 在 确定 模型 参数 个 数 时 ,为 避免 “过 参数 化 ” 
的 出 现 ,将 遵循 规则 :样本 数 和 参数 个 数 的 关系 应 当 
符合 :n 三 3(k +1)(n 为 样本 数 ,k 为 最 终 模型 中 的 
参数 个 数 ) ;相关 系数 随 参数 个 数 的 增加 并 无 显著 
增 大 时 (AR ”<0.02) ,回归 终止 ,此 时 获得 的 模型 为 
最 佳 模型 。 

3) 模 型 的 稳定 性 和 预测 能 力 的 检验 由 “ 留 一 
法 ” 交 又 检验 和 三 重 内 部 检验 完成 。“ 留 一 法 ”交叉 
检验 结果 可 由 CODESSA 软件 直接 获得 。 三 重 内 部 
检验 的 步骤 为 :首先 将 22 个 化 合 物 分 成 3 组 ,将 第 
1,4,7,10 证 等 分 为 A 组 , 2, 5,8 | 等 分 为 B 
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组 , 剩 下 的 为 C 组 。A, B 和 C3 组 进行 两 两 组 合 ， 
获得 新 的 数据 集 A+B, B+C 和 A+C, 利 用 这 3 个 
新 的 数据 集 可 以 构建 3 个 描述 符 个 数 与 最 佳 模 型 相 
同 的 新 的 模型 ,此 时 使 用 的 描述 符 都 选 自 最 佳 模 型 
的 描述 符 ; 最 后 利用 A+B,B+C 和 A+cC 集合 构建 
的 模型 分 别 预 测 C,A 和 B 中 结构 的 活性 ,根据 相关 
系数 判断 其 预测 能 力 。 具 体 方法 参见 Katritzky 等 
(2006, 2008 ) 。 





2 结果 与 分 析 


2.1 羚 酸 与 车 类 驱 避 化 合 物 缔 合 的 能 量 、 距 离 与 
角度 

4 种 羧 酸 引诱 化 合 物 与 驱 避 化 合 物 缔 合 作用 的 
能 量 ,以 及 缔 合 体内 的 相互 作用 距离 与 角度 分 别 列 


于 表 1 ~5。 分析 送 酸 和 驱 避 化 合 物 的 结构 特征 , 同 
时 观察 两 者 相互 作用 的 形式 可 以 发 现 , 凑 酸 与 驱 避 
化 合 物 之 间 主 要 以 氧 键 的 形式 相互 作用 。 由 表 1 数 
据 可 知 ,4 种 羚 酸 - 驱 避 化 合 物 缔 合 作用 的 能 量 平 均 
值 在 18 ~ 20 kJ/mol 范围 之 内 。 由 于 4 种 羧 酸 参 与 
缔 合 作用 的 基 团 相同 ,因此 缔 合 能 量 也 差异 不 大 。 
在 这 4 种 羧 酸 中 ,乙酸 与 驱 避 化 合 物 缔 合 作用 的 平 
均 能 量 略 高 ,可 能 的 原因 是 乙酸 有 较 小 的 分 子 体积 ， 
在 缔 合 过 程 中 ,可 以 更 容易 找到 最 低能 量 构象 。 另 
外 , 丙 酸 、 丁 酸 和 戊 酸 分 别 与 19 号 结构 缔 合 时 ,能 量 
几乎 为 零 ,获得 这 一 结果 可 能 原因 是 19 号 结构 自身 
具有 相对 较 大 的 分 子 体 积 , 在 与 丙 酸 、 丁 酸 和 戊 酸 缔 
合 时 ,由 于 缔 合 位 点 附近 的 位 阻 原因 ,使 得 缔 合 作用 
较 难 发 生 。 观 察 其 他 能 量 低 于 平均 值 的 缔 合 情况 ， 
也 大 多 与 19 号 化 合 物 的 情形 有 类 似 之 处 。 



































表 1 了 驱 避 化 合 物 的 活性 数据 及 其 与 凑 酸 缔 合 的 能 量 


‘Table 1 Repellency of terpenoids and the association energy between terpenoids and carboxylic acids 





引诱 化 合 物 与 驱 避 化 合 物 缔 合 作用 的 能 量 (kJ/mol) 








编号 结构 式 活性 数据 Association energy between attractants and repellents 
No. Structural formula logCRR 乙酸 - 驱 避 物 丙 酸 - 驱 避 物 丁 酸 - 驱 避 物 戊 酸 - 驱 避 物 
Acetic acid-repellent Propionic acid-repellent Butyric acid-repellent Pentanoic acid-repellent 
0 
1 日 1.954 19.6 16.4 20.3 20.2 
OH 
0 
2 CH3 1.908 20.6 17.3 20.1 20.7 
OH 
OH 
3 1.862 15.4 33.8 24.2 33.6 
OH 
0 六 
4 1.860 23.9 24.9 25.2 23; 汉 
0 
5 CoHs 1.857 21.1 19.7 19.6 19.6 
OH 
OH 
6 1.839 2537 18.0 17.9 17.9 
OCOC2Hs 
7 1.809 24.5 24.3 22:7 22.7 
O 
8 0 1.803 21.1 20.3 20.0 20.1 
OC2H5 
9 1.794 18.8 17.3 17.9 17.9 
16.9 16.8 16.8 


OCOCHs 
10 1.784 21.:7 
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续 表 1 Table 1 continued 








引诱 化 合 物 与 驱 避 化 合 物 缔 合 作用 的 能 量 (kJ/mol) 








编号 结构 式 活性 数据 Association energy between attractants and repellents 
No. Structural formula logCRR 乙酸 - 驱 避 物 丙 酸 - 驱 避 物 丁 酸 - 驱 避 物 戊 酸 - 驱 避 物 
Acetic acid-repellent Propionic acid-repellent Butyric acid-repellent Pentanoic acid-repellent 
人 
11 @ CH 1.767 15.9 15.7 15.7 15.7 
12 1.754 16.8 15.8 15.5 15.6 
CH2CH20COC2Hs 
O 
13 人 1.749 17.3 17.1 16.9 17.0 
14 “ 1.724 24.0 23.2 22.9 23.0 
15 1.724 21.0 18.8 18.6 18.6 
CH,CH,OCOCH;, 
16 a T9717 17.6 19.1 19.0 19.0 
17 Se 1.713 15.1 14.3 14.0 14.0 
OCOCH; 
18 Cr 1.712 19.0 20.2 20.1 20.1 
OCOCH, 
19 (CD 1.701 18.9 2.1 0.2 2.0 
2 3 
20 有 1.607 18.1 17.0 16.7 16.8 
0 人 Ca 
O 
21 0 人 1.588 18.1 18.7 18.6 18.6 
O 
22 EE 1.577 17.8 17.8 17.7 17.7 
均值 19.6 18.6 18.2 18.8 
Mean 


表 2 ~5 列 出 了 4 种 其 酸 与 驱 避 化 合 物 缔 合 作 
用 的 距离 和 角度 ,在 缔 合 类 型 中 ,虚线 表示 缔 合 作 
用 ,虚线 前 面 的 基 团 是 由 凑 酸 提供 ,后 面 的 基 团 是 由 
驱 避 化 合 物 提供 。 通 常情 况 下 , 羧 酸 充 当 的 是 氧 键 
缔 合 (0-H 0) 中 氧 的 提供 体 , 驱 避 化 合 物 充当 氧 的 
受 体 。 但 是 ,如 果 驱 避 化 合 物 中 也 存在 羟基 时 ,也 可 
以 作为 氧 的 提供 体 与 凌 酸 进行 缔 合 。 缔 合 距离 的 范 









































围 为 2.0~2.5 A, 而 缔 合 角度 的 在 110° ~ 180° 之 间 
2.2” ”QSAR 模型 

2.2.1 描述 符 的 计算 :经 AMPAC 和 CODESSA 软 
件 计算 ,获得 了 大 量 描述 符 , 单 个 模型 可 以 使 用 的 描 
述 符 数量 达 1 200 多 个 。 更 大 描述 符 集 更 有 利于 得 
选 获得 真正 与 活性 数据 密切 相关 的 描述 符 。 
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表 2 乙酸- 驱 避 物 缔 合 体 的 缔 合 距离 和 角度 


Table 2 ”Association distances and angles of acetic acid-repellent complexes 

















编号 缔 合 类 型 距离 (A) 角度 (°) 编号 缔 合 类 型 距离 (A) 9 度 (°) 
No. Association type Distance Angle No. Association type Distance Angle 
1 -OH …O=C< 2.255 149.053 12 -OH...0=C< 2. 135 170. 826 
-OH … O < 2.400 154.427 
>C=0.… HO- 2.452 162.281 
2 -OH...0=C< 2:157 177.519 13 -OH...0=C< 2.094 169.577 
3 -OH...0< 2.081 146. 406 14 -OH...0< 2.068 153. 194 
>C=0...HO- 2.192 120.370 
-OH...0< 2.099 156.219 15 -OH...0=C< 2.136 168. 906 
5 -OH...0=C< 2.219 144.918 16 -OH...0< 2. 142 176.431 
-OH...0< 2.340 155.913 
>C=0...HO- 2.398 162.921 
6 -OH...0< 2.119 152. 811 17 -OH...0=C< 2.132 173. 802 
>C=0...HO- 2.136 123.407 
7 -OH...0=C< 2;:132 170.088 18 -OH...O=C< 2. 149 165.426 
8 -OH...0=C< 2.133 167. 865 19 -OH...0=C< 2.131 172. 146 
9 -OH...0< 2.095 174.491 20 -OH...0=C< 2.116 173.620 
10 -OH...0=C< 2.106 176. 193 21 -OH...O=C< 2. 103 175.228 
11 -OH...0=C< 2.135 165.414 22 -OH...0< 2:,155 167.071 















































在 缔 合 类 型 中 ,虚线 表示 缔 合 作用 ,虚线 前 面 的 基 团 是 由 凑 酸 提供 ,后 面 的 基 团 是 由 驱 避 化 合 物 提供 。 表 3 ~5 同 。As to the association type， 


dashed line means the association, groups in front of the dashed line are supplied by carboxylic acids and groups behind the dashed line are provided by 











repellents. The same for Tables 3 -5. 


表 3 丙 酸 - 驱 避 物 缔 合 体 的 缔 合 距离 和 角度 


Table 3 Association distances and angles of propionic acid-repellent complexes 

















编号 缔 合 类 型 距离 (A) 角度 (°) 编号 缔 合 类 型 距离 (A) 度 (°) 
No. Association type Distance Angle No. Association type Distance Angle 

1 -OH...0=C< 2.:137 164.477 12 -OH...0=C< 2.136 147.276 

2 -OH...0=C< 2.300 146. 079 13 -OH...0=C< 2.095 169.772 
-OH...0< 2.152 150.449 

3 -OH...0< 2.099 153.506 14 -OH...0=C< 2.067 153.814 

>C=0...HO- 2.158 121.991 >C=0...HO- 2.192 119.954 

4 -OH...0< 2.091 167.614 15 -OH...0=C< 2.097 174.608 

5 -OH...0=C< 2.146 168. 145 16 -OH...0< 2.144 176.375 

6 -OH...0< 2.100 153.498 17 -OH...0=0< 2:127 173.561 

7 -OH...0=C< 2.277 148. 200 18 -OH...0=C< 2.146 165.877 
-OH...0< 2..3357 144.659 

8 -OH...0=C< 2.130 166.915 19 -OH...0=C< 2:127 166.088 

9 -OH...0< 2.104 173.006 20 -OH...0=C< 2.142 175.719 

10 -OH...0=C< 2.247 149.354 21 -OH...0=C< 2.102 175.278 
-OH...0< 2.418 143.740 

11 -OH...0=C< 2.130 165. 824 22 -OH...0< 2.161 166. 292 








表 4 了 丁 酸 - 驱 避 物 缔 合 体 的 缔 合 距离 和 角度 


Table 4 Association distances and angles of butyric acid-repellent complexes 





编号 缔 合 类 型 距离 (A) 角度 (°) 编号 缔 合 类 型 距离 (A) 角度 (°) 
No. Association type Distance Angle No. Association type Distance Angle 

1 -OH...0=C< 2.149 165.616 12 -OH...0=C< 2.117 158.855 

2 -OH...0=C< 2.150 147.562 13 -OH...0=C< 2.095 169. 584 
-OH...0< 2.311 148.684 

3 -OH...0< 2. 161 127.013 14 -OH...0=C< 2.044 153.024 

>C=0...HO- 2.099 155.470 >C=0...HO- 2.178 118.477 

4 -OH...0< 2.090 168. 860 15 -OH...0=C< 2.097 174.432 

5 -OH...0=C< 2.151 169.046 16 -OH...0< 2.143 176. 334 

6 -OH...0< 2.098 154.157 17 -OH...0=0< 2.126 170. 456 

3 -OH...0=C< 2.261 149.922 18 -OH...0=C< 2.147 165. 364 
-OH...0< 2.386 143. 842 

8 -OH...0=C< 2.103 170. 506 19 -OH...0=C< 2.132 170.221 

9 -OH...0< 2.097 173.611 20 -OH...0=C< 2.146 174.921 

10 -OH...0=C< 2.241 150. 568 21 -OH...0=C< 2.102 175. 206 
-OH...0< 2.438 142.991 

11 -OH...0=C< 2.144 169.980 22 -OH...0< 2.158 166.031 























Cn 
cc 
Na 
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表 5 成 酸 - 驱 避 物 缔 合 体 的 缔 合 距离 和 角度 


TableS Association distances and angles of pentanoic acid-repellent complexes 





编号 缔 合 类 型 距离 (A) 角度 (°) 编号 缔 合 类 型 距离 (A) 角度 (°) 
No. Association type Distance Angle No. Association type Distance Angle 

1 -OH...0=C< 2.148 166. 168 12 -OH...0=C< 2.141 145.897 

2 -OH...0=C< 2.147 147.534 13 -OH...0=C< 2.096 169.532 
-OH...0< 2.307 149. 594 

3 -OH...0< 2.096 153.279 14 -OH...0=C< 2.046 153.431 

>C=0...HO- 2:159 121. 102 >C=0...HO- 2.178 119.448 

4 -OH...0< 2.091 168. 840 15 -OH...0=C< 2.100 174.518 

5 -OH...0=C< 2.150 169.449 16 -OH...0< 2.144 176.462 

6 -OH...0< 2.097 154.205 17 -OH...0=0< 2.126 171.442 

7 -OH...0=C< 2.266 149.551 18 -OH...0=C< 2.150 165.412 
-OH...0< 2.383 143. 879 

8 -OH...0=C< 2.136 167.766 19 -OH...0=C< 2.131 169.420 

-OH...0< 2.098 173.572 20 -OH...0=C< 2.143 173.305 

10 -OH...0=C< 2.241 149.744 21 -OH...0=C< 2.102 175.235 
-OH...0< 2.437 143.266 

11 -OH...0=C< 2.135 167.632 22 -OH...0< 2.157 165. 883 


2.2.2 ”QSAR 模型 的 构建 :利用 启发 式 线性 回归 方 
法 筛选 并 构建 了 4 个 最 佳 线性 模型 ,乙酸 - 驱 避 化 合 
物 模型 再 酸 - 驱 避 化 合 物 模 型 、 丁 酸 驱 避 化 合 物 模 
型 和 成 酸 驱 避 化 合 物 模 型 ,模型 的 详细 数据 列 于 
表 6。 

2.2.3 模型 的 检验 :利用 “ 留 一 法 ”交叉 检验 和 三 
重 内 部 交叉 检验 对 获得 的 4 个 模型 进行 检验 。 结 果 
列 于 表 7。 对 比分 析 每 个 模型 各 自 的 相关 系数 平方 
R 交叉 检验 相关 系数 平方 Rs. 训练 集 相关 系数 平 
方 R(tt) 平 均值 和 测试 集 相关 系数 平方 R (pred) 
平均 值 可 知 , 获 得 的 4 个 模型 具有 和 较 好 的 稳定 性 和 
预测 能 力 。 对 于 乙酸 - 驱 避 化 合 物 模 型 ,虽然 模型 只 
含有 3 个 描述 符 , 但 是 模型 检验 的 结果 依然 理想 。 
图 2 所 示 的 是 利用 表 6 中 的 4 个 模型 获得 的 计算 值 
与 实验 测试 值 的 对 比 , 可 以 看 到 ,两 者 之 间 有 很 好 的 
线性 关系 ,相关 系数 平方 都 大 于 0.9。 综 上 所 述 , 获 
得 的 4 个 模型 结果 理想 。 





























3 讨论 


本 研究 利用 CAUSSIAN 软件 计算 了 4 种 羧 酸 与 
驱 避 化 合 物 分 子 间 的 相互 作用 ( 缔 合 作用 ) ,并 借助 
CODESSA 软件 分 析 了 这 种 缔 合 作用 对 驱 避 活性 的 
影响 。 人 研究 获得 了 4 种 凑 酸 与 莫 类 驱 避 化 合 物 相 互 
缔 合 的 缔 合 能 量 、 距 离 和 角度 ,并 获得 了 4 个 最 优 定 
量 构 效 关系 模型 。 从 模型 中 所 包含 的 参数 分 析 得 
出 ,引诱 物 - 驱 避 物 缔 合 作用 对 驱 避 活性 具有 重要 


恒 / 
尿 Z 啊 3 











一 般 而 言 ,分 子 间 氢 键 的 键 能 在 10 ~ 65 kJ/mol 
之 间 , 氢 键 的 距离 根据 构成 氢 键 原子 的 不 同 而 不 同 ， 
氧 键 的 键 角 应 大 于 90° ,接近 180°( 杨 华 等 , 2003 ) 。 
本 研究 获得 的 法 酸 与 驱 避 化 合 物 分 子 缔 合 作用 的 能 
量 平均 值 在 18 ~20 kJ/mol 范围 之 内 ,而 缔 合 作用 
的 距离 范围 为 2.0 ~2.5A, 缔 合 角 度 的 范围 为 
110° ~ 180°; 参 与 缔 合 的 基 团 类 型 主要 为 羟基 恢 基 
等 含 氧 基 团 , 综 上 特征 推断 : 送 酸 与 驱 避 化 合 物 的 缔 
合作 用 属于 中 等 强度 的 氢 键 作用 。 

QSAR 模型 中 的 描述 符 通常 含有 一 些 有 价值 的 
言 息 ,从 而 有 助 于 解释 驱 避 作用 机 理 。 以 下 的 描述 
符 释 义 和 人 性 质 分 析 主 要 参考 自 CODESSA Pro 的 使 
用 参考 手册 (Katritzky et al., 2002 ) 。 

在 乙酸 - 驱 避 化 合 物 模型 中 ,第 1 个 描述 符 :Min 
nucleoph. react. index for a C atom ,为 驱 避 化 合 物 分 
子 描述 符 。 该 描述 符 与 最 高 占用 分 子 轨道 相关 ,是 
琉 电 子 作用 的 指示 描述 符 (Katritzky et al.，2008 ) 。 
它 对 驱 避 化 合 物 的 电子 性 能 可 能 存在 影响 ,改变 这 
一 参数 的 数值 ,可 能 会 改变 驱 避 化 合 物 在 分 子 间 相 
互 作用 的 贡献 度 , 当 驱 避 化 合 物 与 引诱 物 分 子 相互 
作用 时 , 它 可 以 指示 驱 避 化 合 物 接受 孤 对 电子 的 能 
力 。 第 2 个 描述 符 是 与 电子 间 排 斥 能 有 关 , 是 一 个 
量子 化 学 描述 符 。 第 3 个 描述 符 又 是 一 个 与 分 子 轨 
道 相 关 的 描述 符 , 来 自 于 缔 合 体 。 它 意 示 了 最 高 占 
用 分 子 轨 道 与 最 低 未 占用 分 子 轨道 的 能 量 差 ,是 缔 
合体 结构 稳定 性 的 一 个 指标 。 

在 丙 酸 - 驱 避 化 合 物 模型 中 ,第 1 个 描述 符 和 人生 
2 个 描述 符 与 乙酸 - 驱 避 化 合 物 一 样 。 第 3 个 描述 
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表 6 定量 构 效 关系 模型 
Table6 QSAR models 


描述 符 编号 系数 系数 误差 -检验 描述 符 
No. Coefficient Coefficient error t-Test Descriptors 





乙酸 - 驱 避 物 模型 


Acetic acid-repellent model 


R? =0.9065, R?, =0. 8684, RMSE =0.0292, s? =0.0010, F =58.15, n=22, logCRR|, =2.90 -175.04A -0.78A, -0.04A; 


Ao 2.90 0.13 21. 85 Intercept 
Ai -175.04 27. 19 -6.44 mol-Monomer-Min nucleoph. react index for a C atoma 
A， -0.78 0.04 -6.61 mol-Monomer-Min e-e repulsion for a H atom 
As -0.04 0.01 -3.04 mol-TerAA-HOMO-LUMO energy gap? 
丙 酸 驱 避 物 模型 


Propionic acid-repellent model 


R? =0.9530, R2, =0.9264, RMSE =0.0207, s* =0.0006, F =86.16, n=22, logCRR, =2.78 -286.79B, -0.26B, +0.003B; -0.32B， 


Bo 2.78 0.08 37.03 Intercept 

Bi —286.79 21.87 —13.11 mol-Monomer-Min nucleoph. react. index for a C atom 
B, 一 0.20 0.03 —8.30 mol-Monomer-Min e-e repulsion for a H atom 

Bs 0.003 0.0008 —5.84 mol-TerPrA-ZX Shadow° 

Ba -0.32 0.12 —2.75 f-TerPrA-ESP-FHASA Fractional HASA (HASA/TMSA)! 


丁 酸 - 驱 避 物 模 型 


Butyric acid-repellent model 


R? =0.9285, R?, =0. 8802, RMSE =0.0255, s? =0.0008, F =55.16, n=22, logCRRs =0.78 -207.99C -0.28C, -0.002C; -0.15C， 


Co 0.78 0.86 0.91 Intercept 

Ci 一 207. 99 27553 二 9 和 mol-Monomer-Min nucleoph. react. index for a C atom 
C» -0.28 0.04 —7.60 mol-Monomer-Min e-e repulsion for a H atom 

Cs —0.002 0.0009 —2.39 mol-TerBA-ZX Shadow°® 

Ca 0.15 0.07 2.23 mol-TerBA-Max total interaction for a C-H bond 


戊 酸 - 驱 避 物 模 型 


Pentanoic acid-repellent model 


R? =0.9509, R?, =0.9213, RMSE =0.0211, s =0.0006, F =82.37, n=22, logCRRs =2.24 +2.87D -174.63D, -0.83D; +21.72D;, 


Co 2.24 0.13 17.55 Intercept 

Ci 2.87 0.27 10.71 mol-Monomer-FHBCA Fractional HBCA (HBCA/TMSA) 
C» —174.63 19.77 —8.83 mol-Monomer-Min nucleoph. react. index for a C atom 
Cs —0.83 0.18 —4.58 mol-Monomer-ZX Shadow / ZX Rectangle 

Ca 21.72 S527 4.12 mol-TerPeA-Moment of inertia Bf 


amol-Monomer- 是 驱 避 化 合 物 分 子 描述 符 的 前 级 mol-Monomer- is a prefix which means descriptor calculated from repellent itself; bmol-TerAA- 是 乙 
酸 - 驱 避 物 缔 合 体 分 子 描述 符 的 前 级 mol-TerAA- is a prefix which means descriptor calculated from acetic acid-repellent complex;“mol-TerPrA- 是 丙 
酸 - 驱 避 物 缔 合 体 分 子 描述 符 的 前 级 mol-TerPrA- is a prefix which means descriptor calculated from propionic acid-repellent complex; df-TerPrA- 是 丙 
酸 - 驱 避 物 缔 合 体 局 部 描述 符 的 前 缀 f-TerPrA- is a prefix which means descriptor calculated from fragment in the propionic acid-repellent complex; 
*mol-TerBA- 是 丁 酸 - 驱 避 物 缔 合 体 分 子 描述 符 的 前 缀 mol-TerBA- is a prefix which means descriptor calculated from butyric acid-repellent complex; 
fmol-TerPeA- 是 成 酸 - 驱 避 物 缔 合 体 分 子 描述 符 的 前 缀 mol-TerPeA- is a prefix which means descriptor caleulated from pentanoic acid-repellent 
complex. 及 : 相关 系数 平方 Squared correlation coefficient ; R2 : 交叉 检验 相关 系数 平方 Squared correlation coefficient of cross-validation ;RMSE : 
均 方 根 误差 Root-mean-square error; ?2 : 方差 Variance; 下 : Fisher 氏 检 验 Fisher' s test; n: 样本 数 Sample number. 

















符 :ZX Shadow ,是 一 个 几何 描述 符 , 其 含义 是 丙 酸 - 
驱 避 化 合 物 缔 合 体 的 结构 大 小 和 几何 形状 。 第 4 个 
描述 符 的 含义 为 : 缔 合 体 中 参与 氢 键 缔 合 的 所 原子 
电荷 表面 积 占 整 个 缔 合 区 域 表面 积 的 比例 。 它 是 一 
个 电荷 加 权 部 分 表面 积 描 述 符 , 与 表面 电子 分 布 和 
分 子 间 相互 作用 密切 相关 。 

在 丁 酸 - 驱 避 化 合 物 模型 中 ,前 3 个 描述 符 都 和 


丙 酸 - 驱 避 化 合 物 模 型 的 一 致 。 第 4 个 描述 符 是 缔 
合体 碳 氢 键 相互 作用 的 最 大 数量 。 这 一 描述 符 可 能 
与 位 阻 作用 对 缔 合 体形 成 和 稳定 性 的 影响 相关 。 

在 戊 酸 - 驱 避 化 合 物 模型 中 ,描述 符 的 组 成 和 排 
序 略 有 变化 。 第 一 个 描述 符 的 含义 与 丙 酸 第 4 个 描 
述 符 的 含义 类 似 , 意 为 : 氢 键 电 稚 表面 积 占 整 个 缔 合 
区 域 表 面积 的 比例 。 也 是 一 个 电荷 加 权 部 分 表面 积 
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logCRR 的 实验 值 
Experimental value of logCRR 
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logCRR 的 计算 值 Caleulated value of logCRR 
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Experimental value of logCRR ~ 
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logCRR 的 计算 值 Calculated value of logCRR 
图 2 模型 计算 的 校正 驱 避 率 对 数值 与 实验 值 的 线性 关系 

Fig.2 Linear relationship between calculated logCRR and experimental logCRR 

A: 乙酸 - 驱 避 化 合 物 模 型 Acetic acid-repellent model; B: 丙 酸 - 驱 避 化 合 物 模 型 Propionic acid-repellent model; C: 丁 酸 - 驱 避 化 合 物 模型 Butyric 
acid-repellent model; D: 戊 酸 - 驱 避 化 合 物 模 型 Pentanoic acid-repellent model. 
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表 7 模型 的 检验 描述 符 ,其 含义 为 驱 避 化 合 物 分 子 中 参与 缔 合 的 基 
‘Table 7 Validation of QSAR models 团 表面 电子 分 布 情况 。 第 2 个 和 第 3 个 描述 符 的 含 
与 前 几 个 模型 相同 。 第 4 个 描述 符 为 缔 合 体 的 几 
AB 085 CC 0.9573 何 参 数 描述 符 ,表示 缔 合 体 分 子 质量 的 分 布 情况 ,由 
乙酸 - 驱 避 物 AC 0.9554  B 0.8769 于 缔 合 体 是 由 引诱 化 合 物 和 驱 避 化 合 物 两 种 分 子 组 
9 A094 成 ,因此 , 缔 合 体 的 质量 分 布 情况 意 示 两 者 的 相对 位 
Wente Qs 置 ,是 缔 合 方式 和 缔 合 体 几 何 构 型 的 一 个 重要 指示 

AB 0. 8394 C 0.9351 参数 。 

a 此 次 研究 与 之 前 的 研究 相 比 较 ( 雇 圣 良 等 ， 
cd odalient 均值 oon 时 2012 ;Song et al., 2013) ,得 到 了 一 些 相 同 的 结论 , 即 
Mean | I 缔 合 体 表面 的 电子 分 布 和 特定 的 几何 构象 对 于 引诱 
有 上 I 化 合 物 与 驱 避 化 合 物 缔 合 作用 具有 重要 影响 。 在 上 
gps oam Bc oo A oz 述 的 这 些 描述 符 中 , 既 有 来 自 于 驱 避 化 合 物 分 子 本 
acid-repellent 均值 gy ey 号 的 描述 符 , 也 有 来 自 缔 合体 及 缔 合 区 域 的 描述 符 。 
人 RE 在 得 到 的 4 个 模型 中 , 每 个 模型 都 含有 来 自 于 缔 合 
二 二 各 Ke Os 体 及 缔 合 区 域 的 描述 符 , 这 表明 缔 合 体 的 结构 与 驱 
Pentanoic 0.9509 0.9213 BC 0.9737 A -0.8703 避 活 性 密切 相关 。 因 此 ,引诱 化 合 物 与 驱 避 化 合 物 
oe repenen ea 0.9053 O03 的 缔 合 作用 是 影响 驱 避 作用 的 一 个 重要 因素 ,在 驱 





避 作 用 机 理 研究 中 不 容 忽 视 。 
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综 上 所 述 , 本 文 利 用 理论 计算 的 方法 研究 了 
4 种 人 体 引 诱 化 合 物 与 22 个 项 类 驱 避 化 合 物 的 
分 子 间 相 互 作用 ( 缔 合 作用 ) ,以 及 这 种 相互 作用 对 
了 驱 避 活性 的 影响 。 结 果 显 示 这 种 相互 作用 为 分 子 间 
氧 键 相互 作用 ,并 对 驱 避 活性 具有 影响 。 未 来 ,更 多 
的 研究 有 竺 结合 多 种 实验 手段 (理论 计算 、 分 子 生 
物 学 、 触 角 电 生理 .嗅觉 行为 学 等 ) 来 系统 研究 蚊 
忠 、 人 体 、 驱 避 化 合 物 和 引诱 物 之 间 的 内 在 相互 
作用 。 
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蚊虫 饲养 和 化 合 物 分 析 测 定 方面 的 条 件 。 
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